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Abstract 

We study a modified mean-field approximation for the Ising Model in arbitrary dimension. In- 
stead of taking a "central" spin, or a small "drop" of fluctuating spins coupled to the effective field 
of their nearest neighbors as in the Mean-Field or the Bethe-Peierls- Weiss methods, we take an 
infinite chain of fluctuating spins coupled to the mean field of the rest of the lattice. This results 
in a significative improvement of the Mean-Field approximation with a small extra effort. 



1 



The Hamiltonian of the Ising model [1[ in the presence of an external magnetic field B 
can be written as 

H = -JY1 s,s,-BY,s,. (1) 

<ij> i 

The spins Sj can take the values ±1. The first sum is over pairs of nearest neighboring spins. 
The partition function is given by 

z= ^M-PH), (2) 

Configurations 

where (3 is the inverse of the temperature. From now on we absorb the parameter J into f3. 

As it is well known, ([2]) can be exactly summed up only for dimensions d = 1,2 (the 
solution of the two dimensional case is far from being trivial). For higher dimensions we 
must content ourselves with approximate analytical solutions or computer simulations. 

The mean-field (MF) approximation is a very simple and intuitive self-consistent 

approximation method that predicts some rough features of the Ising model, and of other 
models that undergo phase transitions. The method consists on "freezing" all the lattice 
spins to a fixed value h, except one of them (the "central" spin Sq). With this simplification, 
the magnetization per unit of volume < sq > is given by 



<so>= ^(*o) So, (3) 

S0=±1 

where the probability P{so) of finding the value sq for the central spin is 

^ exp{2d(3hso) 
^ °^ exp{2d(3h)+exp{-2d(3h)' 

for the case B = 0, which is the one we are interested in. 

When < So > is equated to the mean-value h, one obtains the transcendental equation 



h = tanh{2dhf3), (5) 

that yields estimates for the magnetization and the critical temperature of the exact model. 
The MF method predicts a phase transition for Pmf = ^, with /i = for /5 < Pmf and 
h ^ for (3 > Pmf- Thus, MF incorrectly yields a phase transition for d = 1. However, as 
d increases the results improve: for d oo, (3mf fits the exact result (see table [T|). 
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The MF approach can be formulated as a variational approximation as follows. Since the 
Gibbs probability distribution minimizes the Helmholtz free energy (p = E — {E and S 
are the average energy and the entropy respectively), the probability distribution ansatz 

P. = ^"P^"^-^-^ , (6) 
Econfig.exp(aEi Si) 

will produce a free energy (pa that overestimates the exact one. Now, varying (pa with respect 
to a it is easy to obtain that its minimum value is attained when a = tanh.{2da(3), which 
just corresponds to the MF ansatz. Hence, (pMF can also be used to estimate the Helmholtz 
free energy (p. Furthermore, it also makes sense to estimate other thermodynamical ob- 
servables from this approximate free energy (pMF- There also exist MF based approaches, 
that instead of considering just one fluctuating spin, take a small sublattice (it could be a 
"central" spin and its nearest neighbors, for instance) and couple it with the mean field of 
the background provided by the rest of the lattice . These approaches improve 

the estimations of the critical temperature and the modeling of the physical observables. 
However, since in essence all of them deal with finite lattices, the estimations of correlation 
functions are somewhat artificial [by considering an external non-uniform magnetic field and 
using the linear response theorem [sl one can estimate correlation functions from the MF 
approximation] . 

Consider now the following approximation, which we propose as a natural improvement of 
the MF one that also allows for a simple estimation of correlation functions. Instead of taking 
a single fluctuating spin, let us consider a straight chain of (A^ ^ oo) fluctuating spins 
in the (i-dimensional lattice, which are coupled to the constant field M of their remaining 
"frozen" nearest neighbors. It is easy to see that this amounts to dealing with the one- 
dimensional Ising model in an external magnetic field B = 2{d — 1)M. The probability 
distribution Pi for this line-approximation is given by 



exp (/5 Ef + /5 Ef 2(t/ - 1)M s,^ 

Pi - y ' (') 

where the partition function Zi is 

I N N \ 

Zi= exp /5^s,s,+i + /5^2(rf-l)Ms, . (8) 

Config. \ i i / 

Zi can be easily calculated using, for instance, the transfer matrix method [Qj. One has 



Tr(f 



(9) 



with 



exp {P + b) exp (— /9) 
exp (— /9) exp (/3 — b) 



(10) 



Here, we have set b = 2{d — 1)(3M. 

The matrix T can be readily diagonahzed. Its eigenvalues are 

A± = exp{(3) (^cosh(6) ± sinh^b) + exp{-4(3)^ . (11) 

In the thermodynamical limit (A^ oo) the partition function only picks the contribution 
of the greatest eigenvalue. The result is 



N 



exp (/3) cosh(6) + Jsinh\b) + exp(-4/5) 



The self-consistency condition (sj) = M is now readily seen to be 



(12) 



M 



d 



2{d-l)(3N dM 
which leads to the transcendental equation. 



InZ; 



(13) 



M 



sinh(6) 



sinh''(6) + exp (-4/3) 



(14) 



Equation f|T4l) always has M = as solution. A second solution with a lesser free energy 
appears for P < Pc, where Pc is given by 



2{d-l)P,exp{2P,) = l. (15) 

Table [T] shows the results for several dimensions. Besides the modest improvement in 
comparison with ordinary MF theory that this approximation yields for intermediate di- 
mensions, the results for the limiting cases d = 1 and d = oo are exact. As in the ordinary 
MF case, the magnetization M vanishes for P < Pc , while M 7^ for P > Pc- In figure [1] we 
show the magnetization vs. P for several dimensions. The behavior is qualitatively similar 
to the ordinary MF case, although the critical temperatures are improved in the present 
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case. From equations (fT^ and (fT5|) one can obtain the magnetization near the critical point 

as 



M 



6(2/?, + 1) 



(5/3) 



/2 + 0{5(3f'\ 



(16) 



^3/3e-exp(-4/5,) 



with 5/3 = /? — /?c > . The magnetization critical exponent is |, just as in the ordinary MF 



case. 



It is worth noticing that the modified MF or line- approximation presented here can 
also be obtained from a variational principle. To see this we consider the trial probability 
distribution 



that may be understood as follows. The whole d-dimensional lattice is partitioned into 
infinite chains parallel to the x-axis. Spins belonging to different chains do not interact. 9 
and a are variational parameters, that must be tuned to produce the best trial probability 
distribution, i.e., that minimizing the trial Helmholtz free energy (p = E — ^S. Here, E is 
the average of the Hamiltonian calculated by using the trial probability distribution, while 
5 = — < InP > is the entropy associated to that probability distribution. Once again, 
using the transfer matrix formalism, the corresponding partition function is readily seen to 
be 



N 



P= U eMOT. 



SiSi+i + a^Si) 



(17) 




Config. X— axes \ i=l i 



= exp(^^) [ cosh(Q;) + sinh^(a) + exp(— 2^^) 




(18) 



where N^j = A^x (number of x-axes), is the total number of vertices of the lattice. 
The Helmholtz free energy is then given by 




(19) 
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When the gradient of (j) with respect to the parameters {6, a) is equated to 0, it is found 
that there is an absolute minimum for 



e = (3 (20) 
2{d-l)[3 d\nZ 

which corresponds precisely to the modified MF solution discussed before. In this case, the 
Helmholtz free energy per site results to be 



(j)/N^ = {d- 1)I3M^ - ln(A+). (22) 

In figure [2] the free energy per site times (3 is depicted for dimension = 3 in comparison 
with the MF one. It is evident that the MF results are improved by using the modified MF 
approximation that we are presenting here. For other dimensions the results are similar, 
and the modified MF approximation also improves the ordinary MF ones. 

The modified MF method also allows us to estimate the correlation length ^, employing 
well known results for the one dimensional Ising model in an external magnetic field. Using 
the transfer matrix formalism it is easy to see that the connected correlation function < 
aiUj >c=< o-iaj > — < cTj >< (jj > is given exactly by 

< (TiCTj >c= constant x exp — ) , (23) 



where the correlation length ^ is given by 



it) 

Here, we have taken the spins at and aj along the same x-axis of "alive" spins. Otherwise 
there would not be correlation among them, just as in the ordinary MF approach. In 
figure [3] we show the correlation length as a function of (3 for dimension d = 3. We observe 
with interest that ^ exhibits a pronounced peak precisely when f3 reaches the critical value 
predicted by the model. For other dimensions the graphics are similar. As it is well known, 
it is expected that correlation lengths diverge at critical points. So, we find it satisfactory 
that this simple approximation can mimic that highly non-trivial behavior in a neat way. It 
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should be stressed that this is an analytical result of the approximation we have discussed, i.e. 
the modified MF approximation, which has no counterpart within the ordinary MF method, 
since in that case there is no correlation among spins (as we pointed out at the beginning, 
however, there are "tricks" to estimate correlation lengths in ordinary MF theory, that 
rely on putting an external non-uniform magnetic field and employing the linear response 
theorem) . 

Summarizing, we have presented a refinement to the MF approximation for the Ising 
Model that: 1) Improves the estimations of the critical temperature, in comparison with the 
ordinary MF approximation, 2) fits exactly to the results of the one dimensional case (unlike 
normal MF) and the infinite dimensional one (as MF also docs), 3) improves the estimates 
for the free energy and, 4) allows to compute correlation functions that exhibit a behavior 
at the critical temperature that strongly resemble the expected behavior of the underlying 
physical models. The approximation is entirely analytical, and heavily relies in the exact 
study of one-dimensional systems. 

It could be interesting to apply this approach to other models whose one dimensional cases 
be resoluble. It also looks reasonable to think in adapting it to the case of lattice gauge 
theories. In the last case, the two dimensional models could serve as the starting point to 
set up the approximation, since gauge-fixing allows to solve exactly the gauge theories in 
this dimension. These and other interesting issues are under work. 
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FIG. 1: Magnetization M versus (3 for several dimensions. 
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FIG. 3: Correlation length versus /3 for d=3. The dashed line corresponds to Pc- 



d 


/3c MF 


(3c modified MF 


Pc other methods 


1 


0.500 


oo 


CXD (exact result) 


2 


0.250 


0.284 


0.441 (exact result) 


3 


0.167 


0.176 


0.222 (Monte Carlo) 


4 


0.125 


0.129 


0.150 (Monte Carlo) 


oo 


1 


1 


1 


2d 


2d 


2d 



TABLE I: Results for fie compared with ordinary MF and other methods. 
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